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Abstract 

The AMANDA neutrino telescope has recently reported the detection of high-energy neu- 
trinos spatially and temporally coincident with the flaring of the TeV blazar 1ES 1959+650. 
At present, the statistical significance of this observation cannot be reliably assessed, however. 
In this letter, we investigate whether circumstances exist where the source can produce the 
flux implied by the coincident events. We show that if the TeV gamma-ray emission observed 
from 1ES 1959+650 or other nearby TeV blazars is the result of accelerated protons interact- 
ing with nucleons, it is reasonable that AMANDA could detect several events during a flaring 
period. Such rates require that the spectral index of the source be rather high (for instance 
~ 2.8 for 1ES 1959+650) and that the Lorentz factor of the jet be fairly small (r ~ 1). 



1 Introduction 



Blazars, the class of Active Galactic Nuclei (AGN) with collimated jets aligned along the direction 
of observation, have been detected over a wide range of wavelengths ranging from radio to very 
high energy gamma-rays. Currently, six of these objects have been identified as TeV emitting 
blazars ^ HI EJ 03 El . These objects are, cosmologically speaking, fairly local (between 140 and 
600 Mpc) and extremely luminous, particularly during periods of flaring activity (~ 10 45 — 10 49 
erg/s, inferred isotropically). Such flares may extend for hours or days at a time (ZHHUHl- 

Blazars typically display spectra with enhanced emission over two energy ranges: the IR/X- 
ray and MeV/TeV peaks. The lower energy peak is generally agreed to be the product of 
synchrotron radiation from accelerated electrons [HJ 1101 111! I12| . The origin of the higher energy 
peak is not yet agreed upon, however. In leptonic models, inverse Compton scattering of syn- 
chrotron photons (or other ambient photons) by accelerated electrons generates this high energy 
emission. In hadronic models [H3 EU E3 E3 E| , on the other hand, MeV-TeV gamma-rays are 
produced through proton collisions with radiation or gas clouds surrounding the object. These 
collisions generate neutral and charged pions which decay, producing very high energy gamma- 
rays and (yet unobserved) neutrinos. The observation of high energy neutrinos would be a strong 
confirmation of the hadronic blazar model. Although roughly 60 blazars have been observed in 
the MeV-GeV range, only 6 have been detected at TeV energies. 

Of the six known TeV blazars, the Markanians are the most nearby [z = 0.031 and 0.033, 
respectively). 1ES 2344+514 and 1ES 1959+650 are not much further (z = 0.044 and 0.047), 
while PKS 2155-304 and 1H 1426+428 are the most distant (z = 0.12 and 0.13). All of these 
except for PKS 2155-304 are northern hemisphere objects, which is necessary to avoid the atmo- 
spheric neutrino background in an Antarctic neutrino telescope such as AMANDA or IceCube 
(PKS 2155-304 could be a potential source for Antares, however). 

During the spring and summer of 2002, 1ES 1959+650 underwent a period of bright flaring 
activity. In May of that year, flaring was observed in the TeV range by the Whipple |18j and 
HEGRA |19| experiments as well as in the X-ray range by the RXTE experiment. The X-ray 
observations of 1ES 1959+650 then showed a gradual decline throughout the following month. 
During that same month, however, Whipple observed the recurrence of flaring in the TeV. The 
observation of this "orphan" TeV flare (with no X-ray counterpart) is in striking disagreement 
with the predictions of the leptonic blazar model. In leptonic models, flaring at TeV energies 
is generally accompanied by simultaneous flaring in the lower energy peak . This argument 
designates 1ES 1959+650 as a TeV blazar likely to generate gamma-rays hadronically |21j . and 
therefore likely to generate a substantial flux of high energy neutrinos. 

Very recently, the AMANDA collaboration has reported the detection of two neutrinos coin- 
cident with TeV flares seen by Whipple from the direction of 1ES 1959+650 |22|I23|. These events 
were not uncovered in a blind analysis, however, and therefore their statistical significance cannot 
For arguments attempting to reconcile 1ES 1959+650 with leptonic blazar models, see Ref. |20| 
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be reliably estimated. The AMANDA collaboration has adopted the position that a potential 
signal, studied after the data has been unblinded, can no longer be statistically evaluated. 

Although the probability of background events occurring coincidentally with the orphan 
flare of 1ES 1959+650 is quite small (~ 10~ 3 ), we emphasize that it is impossible to treat this 
statistically, as there is no known number of trials to dilute the result by. With this in mind, 
we proceed to explore the implications of these events, in the case that they are the product of 
hadronic interactions in the jet of 1ES 1959+650. 

Our conclusion is that the events observed from the direction of 1ES 1959+650 by the 
AMANDA experiment are consistent with the model of proton-proton collisions in the jet of a 
hadronic blazar with the characteristics of 1ES 1959+650. 



2 Inferring a Neutrino Spectrum From Gamma-Ray Observa- 
tions 

Assuming the TeV emission observed in a blazar's spectrum is generated in the decay of pions 
produced in the collisions of protons accelerated in its jets, high energy neutrinos must also be 
present. Furthermore, the gamma-ray and neutrino fluxes can be related by energy considerations 
[2H: 

E^dE, = K E u ^dE u , (1) 

where K is a factor which depends on whether the pions are generated in pp (K = 1) or p'j 
{K = 4) collisions. 

The spectrum of protons generated via Fermi acceleration in the blazar jets is expected to be 
of the form dN p /dE p cx E~ 2 , with both the gamma-ray and neutrino spectra tracing this slope 
(assuming a slowly varying spectrum for the target photons/protons). Gamma-rays can interact 
in the source, as well as with the Infra-Red Background (IRB) during propagation, however, 
steepening their observed spectrum. The threshold for the process, 7 + 7irb — ► e + e _ , is given 
by: 

2m 2 

E -r > f M 6 m ~10-50 TeV. (2) 

£i RB (l - cos0) 

Here, -Eire ~ 0.01 eV is the energy of a typical IRB photon and 9 is the angle between the 
two photons. The precise behavior of this interaction on the spectrum of high energy photons 
depends on the poorly measured distribution of extragalactic infrared radiation, as well as on 
the distance to the source. 

Within the source itself, gamma-ray cascading can also significantly modify the observed 
spectrum. These processes make it impossible to match the shape of the neutrino spectrum to 
that of the gamma-rays. The total energy between photons and neutrinos can still be related, 
however. Even this only yields a lower limit on the neutrino flux, however, as absorbed photons 
are not accounted for in the left hand side of Eq. ^ 
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2.1 Proton-Proton Collisions 



Following simple relativistic kinematics, for the center-of-mass energy of a proton-proton collision 
to exceed the energy threshold for pion production, the accelerated proton must have an energy 
above: 

(2m v + m^) 2 — 2ml 

E™ m = r ~ r x 1.23 GeV, 3) 

y 2m p 

where T is the Lorentz factor of the jet relative to the observer. The maximum energy to which 
protons are accelerated depends on the kinematics in the jet. The maximum and minimum 
gamma-ray and neutrino energies are related to this quantities by: 

rmax rpm&x 

' 6 ' v 12 ' { ' 

with analogous expressions relating the minimum energies. The factors of 6 and 12 come from, 
on average, three pions being produced and each charged (neutral) pion decaying into four (two) 
particles. 

2.2 Proton-Photon Collisions 

Pions can also be generated in pj collisions via the delta hadron resonance, p"f —* A — ► 7riV. For 
this process to take place, the center-of-mass energy of the interaction must exceed the A-mass, 
1.232 GeV. This correspond to an energy of: 



) x 160 GeV. (5) 



This is clearly a much higher energy cutoff than in the case of proton-proton collisions. To 
generate TeV photons by this mechanism, radiation clouds surrounding the blazar must contain 
photons of ~MeV energies EH HE] • Such clouds are known to exist around at least some 

line-emitting blazars. 

In the case of pj collisions, the maximum and minimum neutrino and gamma-ray energies 
are given by: 

pmax ^ „ v. pmax ^ ,.. ^ 

irnmax ^P ^p^n max & p ^p^ir ■> , . 

^7 = o ' = a > W 



again with analogous relations for minimum energies. < x p ^ n >~ 0.2 is the average fraction of 
the proton's energy which is transfered to the pion. 

3 Estimates of High Energy Neutrino Fluxes 

Consider a TeV blazar with the observed spectrum: 
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The expectation for the gamma-ray spectrum from Fermi accelerated proton interactions is a ~ 2, 
although the observed spectrum is likely to be steepened due to cascading. We will assume that 
the neutrino spectrum more closely mimics the Fermi accelerated proton form: 



dE v 

Using Eq. ^ we can relate these coefficients: 



A V E~ . (8) 



A p-Q+2 

A ~ 7 7 ' min (Q) 

,min ) 

where we have assumed that i? 7i max 3> Ej >m i n . The neutrino spectrum from a hadronic TeV 
blazar with the spectrum of Eq. [7| is thus given by: 



~ A P-2 _ ^^mtn ,, 

dE v ~ v » ~ (a-2)Kln(E ViXDax /E u ' ' 



which for proton-proton collisions is approximately: 

dN v o i 2m ,r 21/ A 7 \ /10\ Q ~ 2 /2.8 - 2\ . , 

— -El~ 2 x 10~ 12 TeV cm" 2 s" 1 x -2- — r — , (11) 

d£7„ " V3.2 x 10- n TeVcm- 2 s-V V E 7 U"2 / 

and for proton-photon collisions is: 

dN » F 2 zlvin-^TpVrm- 2 ,-^^ ^ \/10\ 2 "- 4 / 2.8-2 \/ E 7 , targct \"- 2 

_^~4xl0 TeV cm s x ^ - ^ — —g -j j ^-—^ ( vTW J • 

(12) 

In these expressions, T is the Lorentz factor of the jet. The gamma-ray flux is normalized to that 
of the Crab nebula: dN^/dE^ ps 3.2 x KT 11 TeV cm" 2 s _1 at 1 TeV. 

At the source, two thirds of the neutrinos produced are of muon flavor (the remaining being 
electron neutrinos). After oscillations which occur during propagation, the neutrino flux reaches 
Earth in roughly equal fractions of electron, muon and tau flavors. 

4 Event Rate Estimates 

To detect a (muon) neutrino of GeV-TeV energy in a high energy neutrino telescope, an energetic 
muon must be produced in a charged current interaction near the detector. Once generated, the 
muon loses energy at a rate: 

dE„ 



— H fa -A- BE U: (13) 
dX ^ K J 

where the values of A and B depend on the medium. For ice, A ~ 2.0 x 10~ 3 GeV cm" 1 and 

B ~ 4.2 x 10 -6 cm -1 . The distance a muon of an energy, E^, travels before its energy drops 

below an energy threshold, £^ hr , is then given by |29| : 



A + BE^ 



A + BE^ 



(14) 
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This quantity is called the muon range. The number of neutrino-induced muons observed in a 
high energy neutrino telescope is given by: 



where A e ^ is the effective muon area of the detector, T is the length of time observed (or the 
duration of a flare) and P u ^^ is the neutrino to muon transition probability, given by: 



where 0cc,^n is the charged current neutrino-nucleon cross section, is the muon range and n 
is the number density of nucleons in the target material. 

For a specific experiment, the event rates predicted are determined by sophisticated monte 
carlos which take into account the detailed geometry and other characteristics of the detector. 
For the AMANDA detector, the effective neutrino area for a northern hemisphere source is 
about 1, 100 and 3000 square centimeters for neutrinos of 100 GeV, 1 TeV and 10 TeV energies, 
respectively. This effective neutrino volume function, convolved with the incoming neutrino 
spectrum, produces the observed number of events. 

In figures^ 31 ^ [21 we plot the predicted number of events in AMANDA from hadronic blazar 
flares for proton-proton and proton-photon collisions, respectively. For each curve plotted, we 
normalize the rate by setting the gamma-ray flux integrated above 2 TeV to that of the Crab 
nebula, F^{E 1 > 2 TeV) sa 5.6 x 10 -12 photons cm -2 s" 1 , multiplied by a 90 day exposure (flare 
duration) . 

During its flaring activity in 2002, the blazar 1ES 1959+650 was observed by two Atmo- 
spheric Cerenkov Telescopes (ACTs), HEGRA and Whipple. The HEGRA collaboration reported 
a series of flares occurring between May and September of 2002 with integrated fluxes often near 
or exceeding the flux of the Crab nebula. The brightest of these, which occurred in May, was 
measured to have an integrated flux above 2 TeV of 1.2 x 10~ n photons cm -2 s _1 , about twice 
the Crab flux. The HEGRA collaboration estimated a spectral index of 2.83 ± 0.14 above their 
threshold of 2 TeV during this activity ^H] • 

The Whipple telescope, operated by the VERITAS collaboration, also reported a series of 
flares during this period, the brightest of which being almost four times the Crab flux (above 
600 GeV). The spectral index reported by Whipple, 2.82 ± 0.15 ± 0.3 [HOI, is in good agreement 
with the value estimated by HEGRA. 

Inserting these values for the 2002 flaring of 1ES 1959+650 into Eqs. ^2 an d 1121 we can 
estimate the neutrino flux and event rates. For proton-proton collisions, using a = 2.83 and 
Ay = 2 x 10 -5 , we estimate a minimum of 1.8 events observed by AMANDA for a Lorentz factor 
of r = 1 or 0.3 and 0.06 events for T = 10 and 100, respectively. The number of events may be 
considerably higher if a substantial fraction of the energy in gamma-rays is absorbed or cascaded 
below the energy thresholds of HEGRA and Whipple. 




(15) 



(16) 
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Figure 1: The lower limit on the number of events in AMANDA predicted from a hadronic blazar 
flare which produces TeV gamma-rays via proton-proton collisions as a function of the Lorentz 
factor of the jet relative to the observer. The solid, dotted and dashed lines correspond to gamma- 
ray spectral indexes, a, of 2.2, 2.5 and 2.8, respectively. The results are normalized to a gamma- 
ray flux integrated above 2 TeV equal to that of the Crab nebula, Fj(Ej > 2 TeV) ~ 5.6 x 1CP 12 
photons cm -2 s _1 , multiplied by a 90 day flare duration. E^ x and S™ ax were set to 20 and 100 
TeV, respectively. Note that the rates have not been increased for absorption of photons on the 
infrared background, nor for possible absorption of the flux in the source. 

For the case of proton-photon collisions, the neutrino rates are considerably smaller, i.e. 0.01 
for r = 1 and MeV target photons. The proton-photon collision scenario appears to be beyond 
the reach of current detector technology (unless the source is not transparent or absorption on 
IR background very large), but may be within the sensitivity of next generation experiments, 
such as IceCube. 



5 Discussion and Conclusions 



The recent observation of neutrinos from the direction of the TeV blazar 1ES 1959+650 coincident 
with flares occurring in 2002 raises the question of whether such large neutrino fluxes can be 
generated in these objects. We have found here that such rates are indeed possible for a hadronic 
blazar which i) accelerates protons which produce pions on nucleonic, rather than photonic, 
targets, ii) the Lorentz factor of the jet if fairly small (r ~ 1), and hi) the observed spectral 
index is fairly large. Given the orphan flare observed in 1ES 1959+650, this object is one of the 
most probable candidates for a gamma-ray source of hadronic origin. It is also not difficult to 
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Figure 2: The lower limit on the number of events in AMANDA predicted from a hadronic blazar 
flare which produces TeV gamma-rays via proton-photon collisions as a function of the energy of 
the target radiation divided by the Lorentz factor of the jet with respect to the observer squared. 
The solid, dotted and dashed lines correspond to gamma-ray spectral indexes, a, of 2.2, 2.5 and 
2.8, respectively. The results are normalized to a gamma-ray flux integrated above 2 TeV equal 
to that of the Crab nebula, F 1 (E 1 > 2 TeV) « 5.6 x 10~ 12 photons cm~ 2 s _1 , multiplied by a 90 
day flare duration. £™ ax and £™ ax were set to 20 and 100 TeV, respectively. Note that the rates 
have not been increased for absorption of photons on the infrared background, nor for possible 
absorption of the flux in the source. 

imagine a mild boost factor in its jet. Finally, the spectral index of 1ES 1959+650 observed by 
HEGRA is 2.83, which satisfies the third criterion. 

It is clear that relatively large rates are a consequence of the steeper spectrum for photons 
relative to neutrinos. One may ask whether the column density in the source required to change 
a spectral index of -2 into -2.8 is reasonable. Steeper spectra implies higher luminosity given 
by the left hand side of Eq. 1. Here again the answer is affirmative. Our highest event rates 
correspond to luminosities of 10 -10 48 erg/s, not unreasonably high for a flaring blazar. Our 
calculation also illustrates that where neutrino emission is concerned, not all blazars are alike 
and that environmental conditions in the jet can be the origin of large variations in the neutrino 
flux. 
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